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Abstract 

This supplement illustrates application of adaptive observer design from [2] for systems which are not uniquely identifiable, ft 
also provides an example of adaptive observer design for a magnetic bearings benchmark system [f ] . 



Key words: Adaptive observers, nonlinear parametrization, weakly attracting sets 



> 

O 
(N 

O 

o 
m 



% 



1 Application of the design to non-identifiable 

systems 

Example 1.1 Consider the following system 



x = Ax + B6iA + g(y), 
2/ = C^x 



(1) 



where 





liy) 



and let ile = [1,3], 51a = [0.5,3], to = 0. In this par- 
ticular case one can, in principle swap and A (and rt\ 
with fig respectively). It is clear that Assumption 3.1 
is satisfied for this system and y{t), as a function of t, 
is bounded for all t > Iq. Moreover, the function ip : 
f{X) = X and the function g arc both continuous and 
satisfy Assumption 3.2. 
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Let us now move to Assumption 4.1. Condition Al is 
satisfied for the given parametrization. Notice, however, 
that if the domain fix contains then the condition 
would not hold. In this situation if docs not belong to 
the domain fig then swapping the definition of A and 9 
(and fix with fig) resolves the issue. Consider condition 
A2. Sets £ and £q, as follows from their definition, coin- 
cide for this system (matrix G is orthogonal to B), and 
are defined as follows: 

£:o(A, 9) ^ £{X, 9) = {(A', 9'), A', 9' e Rj A6' - X'9' = 0}. 

We need to check if there is a function (3 such that 

|A6i - X'9'\ > I3{dist ((A', 9'),£{X, 9))). (2) 

In view of Remark 9, it is enough to show that for any 
{X,9) e fl\X fig there is a function /3 (possibly dependent 
on A, 9) such that (2) holds. And in fact it is enough to 
showthat it holdsforall A', 0' from the domain 17^ xSlg to 

which the estimates A, 9 may belong for t > tg. It is easy 
to see (from the proof of the theorem) that Assumptions 
3.1, 3.2, and Al from Assumption 4. ensure that both 
estimates A, 9 are bounded for t > Iq. Moreover, for the 
given set of initial conditions for 9, one can estimate 
fl*x X fig a-priori. In this example we set fl*^ x fig — 
[0.5,3] X [-10,10]. 
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Let, for the sake of certainty, set A = 2, 6 = 1.5. In 
this case the set £ consists of two isolated curves. For 
any given point (Ae,6'e) from Q,\ x $7^ one can derive 
two projections to the set S. These projections wih be 
on the hnes intersecting each branch orthogonally and 
passing through (Ae, 9e)- Let d^ be the minimal distance 
from (Ae, 9fJ) to these projections. It is clear that for this 
distance de there will be a compact set V^ of points from 
Q,\ X Q,*g from which the distance to £ is exactly d,, . Given 
that De is compact and that \\0 — \eOe\ is continuous 
one can always define 

§: ^(4) = niin \\e - XeOel 

The function f3 can now be defined as a function from 
K-oa satisfying: /3(4) < i?(4)- 

If finding the function d analytically is complicated then 
it can be estimated numerically. To illustrate the proce- 
dure we randomly sampled n^ x Wg and derived the val- 
ues of jA^ — Ae6'e|, de for each pair (Ae, 6*6). The values of 
\\6 — \eOe\ then were plotted against de (please see the 
left panel in 1). The lower envelope of the resulting set of 
points would then be an estimate of d. One can observe 
that in this particular example this envelope is already a 
strictly monotone function. This clarifies condition A2. 

Finally, we simulated the system together with the ob- 
server (with 7g = 1, 7 = 0.1, I. = (—2, —1)^) and the 
results are summarized in Fig 1. 

Notice that the vector of the estimates, (A, 6), may con- 
verge to points outside of the domain VL\ x J7e. It is, how- 
ever, guaranteed that \{t) G Q,\ for all t > tg. Moreover, 
if unmodeled dynamics is present, they may converge 
(depending on initial conditions) to neighborhoods of el- 
ements from £ that are quite far from rix x fig 

Example 1.2 Consider now system 



X = Ax + B6' 



;(j/,A), 



g(2/,A) = 




where A, B, and C, ilg, fix, A, 6 are defined as in (1). It 
is clear that it satisfies Assumptions 3.1, 3.2, and con- 
dition Al of Assumption 4.1. Let us check condition A2 
of Assumption 4.1. First, we notice that sets £ and £o 
are different in this example. Indeed the set £o consists 
of the single element, (A, 6), whereas the set £ is 



£{x,e)^{{x',e'), A', 



0} 



Similarly to the argument presented for the previous ex- 
ample, one can conclude that the function /? (possibly 
dependent on A, 9) exists, and can be estimated numeri- 
cally (see Fig. 2) We simulated the system together with 
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Figure 1. Top panel: scatter plot of 10000 pairs 
(|A6' — Ae6'e|,de) when {Xe,6e) are sampled from the set 
[—10, 10] X [0.5, 3]; red curve shows a boundary of this set 
taken as an estimate of the function j3. Central panel: red 
circles show the values to which the estimates A, 9 converge, 
green circles are the initial conditions, Ao, 6o, and the cor- 
responding branch of £ is shown as the blue curve. Bottom 
panel: variable X2 and £2 as functions of t. 



I 

CD 

I 

+ 



^W"' '' 


y^ 




y^ 










§*"•-' tw ' 


>^ 


fe\ * 


j/^ 




^/^ 


S^;t;* 


<:i--p 


^;r / 


^i^-^'^'^-l'" ■■■ .y^ 


" 


S^^'S'.^y?..' .'■ . .: ■ -■ ' "j^ 






- 




- 



0.4 0.8 1.2 1.6 d^ 



0.5 



25 



15 




1.5 



2.5 




10 



20 



30 



Figure 2. Top panel: estimation of the function /? (obtained 
from 40000 pairs of points {\e + e^ - 8^ - e^^\,dj}; (Ae, 9^) 
are sampled from the set [—10, 10] x [0.5, 3]). Center panel: 
results of the estimation. Bottom panel: X2 (blue) and X2 
(red) as functions of t. 



in Fig. 2. It is evident from the figure that the state com- 
ponent X2 is reconstructed with an error, as expected 
(the Theorem guarantees reconstruction of state only if 
£ and Eq coincide) . 



2 State and parameter estimation for the mag- 
netic beam benchmark system [1] 

The system dynamics is governed by the following equa- 
tions [1] (see also Table 1, example 2): 



x= ( ° ^ )x + ^(i,y(t),A)0 + g(i), y = xi, (4) 
\0 / 



where 



*(i,2/W,A)- r)'^(i,2/(t),A); g(t) = 



(0.5t) 



1000 



(/)(i,2/(t),A) - '^j \(a(y(t)L+g) + \by' (a(-y(t)L+g)+Xbr ) ' 

variables qi (t) , (72 (t) are defined as 
qi = -§hqi + -^sat(ui) 



<72 = -§hq2 + -^sat(u2) 



sat(M) — sign(M) minjy^, |m|}. 



and, h, the linear term of the otherwise nonlinear fiux- 
current dependence, is modeled for simplicity as con- 
stant, h = 10^. Parameters a, b, c, R, N, L, J, Vg are de- 
fined as in [1], and stabilizing control inputs ui,U2 are 
chosen as: 

2x1 -I- 2x2 > ^ 
ui = 



U2 = N{^ - 100(92 
2x1 + 2x2 < ^ 

U2=0 

ui = 7V(^ - 100(91 



2X1+2X2.1 ^ „2 
cL 



qf)) 



(5) 



23:i+2a2./ 
cL 



+ 9I)) 



-20(xi — Xl) + X2 

-lOO(xi-xi)-(-0(i,2/(t),l) 



(6) 



The design is a combination of the high-gain observer (6) 
with the switching controller (5) aimed at ensuring that 
solutions of the unperturbed system obey linear second- 
order equation with asymptotically stable dynamics (cf. 

[ID. 



the observer (7^ = 1, 7 = 0.05, £ = (—2, —1)"^) for vari- 
ous initial conditions in Qx x fig. The results are shown 



Parameters 0, A account for deviations of the true values 
of the magnetic constants a, b, c, and the air gap length, 




Figure 3. Estimates 9, A of S, A and true values of 6, A (black 
solid lines) 

(7, from their nominal values used in (5), (6). In particu- 
lar, we suppose that the true values of 0, A are unknown 
and belong to the interval [0.8, 1.2]. 

Let us now construct an observer for asymptotic re- 
construction of 9, A. Equations (4) are of the type (7), 
(71), and hence the observer candidate is defined by 
(72), (28)- (30). In particular, let B = (1, 1)'^ and define 
M(i,[A,y]), Cas: 



References 

[1] Z. Lin and C. Knospe. A saturated high-gain contol for 
a benchmark experiment. In Proceedings of the American 
Control Conference, pages 2644-2648, 2000. 

[2] I. Tyukin, E. Steur, H. Nijmeijer, and C. van Leeuwen. 
Adaptive observers and parameter estimation for a class of 
systems nonlinear in the parameters. Automatica, 2013. 
(preprint available at http://arxiv.org/abs/0903.2361). 



M 






Vo-1 

' /O 0\ ~ / -2 

+ Bm2iit,[X,y])0 + g{t) 
0' = -7e((l,O)C-2/)m2i(t,[A,y]) 
A = 0.8 + 0.4(si,i + l)/2, 





<j>{t,y{t),X) 

((i,o)C-y) 



(7) 



where the variable si_i is defined as in (79) with 
e = (since no measurement noise and unmodeled 
dynamics are present). One can check that the vec- 
tor {■m2i{t,[X' ,y]),Ri{t,X' ,X,9)) (see the remark after 
Theorem 13), where 

Ri{t,X',X,e) = 

/o Ito e"'*""' i'^(^' 2/(t), s(e, X',X))ddTdC, 
s(C,A',A) = A'e+(l-e)A, 



is A'-UPE for A' G [0.8,1.2], 6 G [0.8,1.2], and hence 
the system is uniquely identifiable, and assumptions of 
Theorem 13 hold. System (4) together with observer (7) 
was simulated for various initial conditions with jg — 
10^, 7 = 150, and samples of typical trajectories are 
shown in Fig. 3. 



